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ABSTRACT 

We study the environments of 6 radio galaxies at 2.2 < z < 2.6 using wide-field near-infrared 
images. We use colour cuts to identify galaxies in this redshift range, and find that three of 
the radio galaxies are surrounded by significant surface overdensities of such galaxies. The 
excess galaxies that comprise these overdensities are strongly clustered, suggesting they are 
physically associated. The colour distribution of the galaxies responsible for the overdensity 
are consistent with those of galaxies that lie within a narrow redshift range at z ~ 2.4. Thus the 
excess galaxies are consistent with being companions of the radio galaxies. The overdensities 
have estimated masses in excess of 1O 14 M , and are dense enough to collapse into virizalised 
structures by the present day: these structures may evolve into groups or clusters of galaxies. 

A flux-limited sample of protocluster galaxies with K < 20.6 mag is derived by statisti- 
cally subtracting the fore- and background galaxies. The colour distribution of the protoclus- 
ter galaxies is bimodal, consisting of a dominant blue sequence, comprising 77 ±10% of the 
galaxies, and a poorly populated red sequence. The blue protocluster galaxies have similar 
colours to local star-forming irregular galaxies (U — Vab ~ 0.6), suggesting most protocluster 
galaxies are still forming stars at the observed epoch. The blue colours and lack of a dominant 
protocluster red sequence implies that these cluster galaxies form the bulk of their stars at 
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1 INTRODUCTION 

In the local Unive rse galaxy clu sters are dominated by red 
early-type galaxies jDresslejri980|) that sit on a tight linear se- 
quence in colour-magnitude s pace known as the red sequence 
(Visvanathan & Sandagdll977t) . These galaxies contain passively 
evolving stellar populations, and have little on-going star for- 
mation. The i ntrinsic scatter of g alaxies about the re d sequence 
is small (e.g.. iBower et alj|l992l ; Istanford et alj|l998h . and both 
the sc atter and slope of the sequence does not e volve with red- 
shift jBlakeslee etafl 120031 : iGladders et ail 1 19981) . These results 
suggest that the stars in cluster elliptical galaxies formed in a 
brief but vigorous epoch of star formation. Current estimates 
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place this epoch of intense star formation at z = 2 — 3, which 
we can th us assume to be an important time in the formation 
of clusters l van Dokkum & van der Marelll2007l : iKurk et alj[20Qg| : 
IPapovich et alj20ld : lTanaka et al.ll201Ch . 

A protocluster is a large-scale overdense region in the early 
Universe, whose properties are consistent with it being a progenitor 
of a local galaxy cluster. Such objects are observed at a time when 
their member galaxies are not yet virialized within single cluster- 
sized dark matter halos. To qualify as a protocluster, a given over- 
dense region must have a mass comparable to that of local clusters 
of galaxies (M > 10 M^), and be dense enough to collapse before 
z = 0. 

An efficient technique for finding the progenitors of clusters 
involves targeting high redshift radio galaxies (HzRGs). Galaxy 
formation is expected to be more efficient in high density en- 
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vironme nts than in the field , due to the abundance of surround- 
ing gas dGunn & Gotdl 1972h . Therefore one may expect HzRGs, 
which are among the most massive galaxies in the ea rly Universe 
dRocca-Volmerange et alj|2004l ; ISevmour et alj|2007h . to pinpoint 
the location of the highest density regions, such as protoclusters. 
Searches for galaxy overdensities near these objects have resulted 
in detecti ons of large-scale structures and protoclusters up to z = 
5.2 (e.g., iLe Fevre et alj 1 19961: iPentericci et al .1 I l2000l: iKurk et al 



200(1 [20041; iBest et alj|2003l : l0verzier et al.ll2006l ; IVenemans eta! 
2007lJoverzier et alfeOO^T 

We have performed an infrared survey of the fields of 8 pow- 
erful radio galaxies with redshifts in the range 1.7 ^ z ^ 2.6 us- 
ing the wide-field near-infrared imager HAWK-I on the VLT. The 
aim of the survey is to identify protoclusters and study the member 
galaxies and their evolution over this redshift range. The HAWK- 
I data on the two lowest redshi ft targets MRC 1017 — 220 and 
MRC 0156 - 252 are presented in lGalametz et alj d2010h . 

In this paper we present HAWK-I data for the 6 highest red- 
shift targets, ranging from 2.2 ^ z ^ 2.6. In Section[2]we describe 
the observations and data reduction. Sections [3]-|5]lay out the evi- 
dence that 3 out of 6 of the radio galaxies lie within galaxy proto- 
clusters. In Section|6]we discuss the properties of the overdensities, 
and the colours of the galaxies within the 3 protocluster candidates. 

Fluxes are calibrated using the Vega magnitude scale unless 
noted otherwise. A flat ACDM cosmology is assumed throughout, 
withfl M =0.3, £2 A = 0.7 and H =70kms _1 Mpc~' . The distance 
scale at z = 2.4 is 8.139 kpc/", or ~ 1.7 co-moving Mpc/'. Distances 
are given in co-moving units, unless stated otherwise. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Target selection and observations 

Six HzRGs with redshift s between 2.28 and 2.55 w ere selected 
from the compendium of iMilev & De Breuckl J2008h . The targets 
were chosen based on their distributions in right ascension, decli- 
nation, and redshift, and their bright radio luminosity (L500MHZ > 
10 28 ' 5 WHz~'), and not on any prior information about surround- 
ing galaxy overdensities. Co-ordinates of the selected targets and 
the control field are given in TableQ] 

The 6 HzRG fields and a blank control field were observed 
in service mode using the High Acu ity Wide field K-band Imager 
(HAWK-I; iKissler-Patig e7ail2008l) on the ESO Very Large Tele- 
scope (VLT) UT4-YEPUN telescope in Paranal, Chile during the 
period April-September 2008. HAWK-I is a near-infrared camera 
comprising of four Hawaii-2 2048x2048 pixel detectors separated 
by a gap of ~15". The camera spans 7.5 x 7.5 arcmin with a pixel 
scale of 0. 106 arcsec per pixel. 

The telescope pointing was optimised so that the HzRG were 
placed close to the centre of the HAWK-I field of view, but at least 
1 arcmin away from the chip gaps, and bright stars did not fall 
within the field of view. Each target was observed through the J, H 
and Ks filters and total exposure times for all targets are provided 
in Table.Q] The telescope dithered every 2 mins which resulted in a 
~1 arcmin- wide cross-shaped region on the final mosaic where the 
image depth is shallower. 

During the observing period the anti-reflection coating of the 
Dewar window of HAWK-I was damaged resulting in small cross- 
shaped patterns on the exposures. These patterns, caused by the 
spider of the secondary mirror, rotated as the telescope tracked the 
targets as HAWK-I is situated on the Nasmyth focus. These pat- 
terns were not adequately removed by the background subtraction 



step unless the telescope had not moved far between successive 
exposures. Thus the integration time between dithers was reduced 
from 2 mins to 1 min for all images obtained after May 2008. 



2.2 Data reduction 

The data were reduced using the ESO/MVM dVandamef2 004) data 
reduction pipeline optimised for the reduction of our HAWK-I data. 
The usual near-infrared reduction steps were taken, including dark 
subtraction, flat-field removal, harmonising the gain of the four de- 
tectors, removal of fringing, sky-subtraction, creation of bad pixel 
and weight maps, calculating the relative astrometry between the 
chips and the absolute astro metry. 

USNO-B1 catalogues jMonet et alj|2003h were used to cal- 
culate the first-guess astrometric solutions for the Ks images. A 
catalogue was then compiled from the objects detected in the Ks 
images and used to calculate the astrometric solution of the J and 
H images. The accuracy of the relative astrometry between the 3 
images of each target is typically within 0. 1 pixel, but the absolute 
astrometry is limited by the USNO-B1 catalogue which has an ac- 
curacy of 0.2 arcsec, equivalent to 2 HAWK-I pixels. 

After the sky background was removed from each exposure by 
MVM in a two step process, low level large-scale variations were 
still seen across each quadrant and most notably between the quad- 
rants of the reduced images. Therefore the sky background of the 
final images was measured and subtracted using a local background 
estimator with SExtractor (BACKPHOTO.TYPE LOCAL). 

The data were flux calibrated with 2MASS catalogues, using 
13 — 15.5 mag stars within the target fields. The calibration was 
checked using 13 — 14 mag stars within standard star fields taken 
within 2 hours of the observations. Typically the zero-points de- 
termined from both methods agreed within the uncertainties (~ 
0.05 mag). 

The end products of this reduction process are a science image 
containing the reduced data of the target, and an effective exposure- 
time image, which is an exposure-time map that has been nor- 
malised to account for differences in sensitivity between the four 
detectors of HAWK-I and the chip gaps. 

The 3 colour images of a target were convolved to match the 
lowest resolution image using the IRAF package PSFMATCH. The 
convolution kernels were optimised so that the stellar growth curves 
(created by median combining at least 20 bright and unsaturated 
stars in each image) converged to within 3% at and beyond a 1" 
radius. The 5o" image depths given in TableQ] were measured by 
placing 2" diameter apertures at multiple random positions. 

2.3 Masking of cross-talk, bright objects, and regions with 
low exposure time 

Each of the four HAWK-I detectors contain 32 amplifiers. Cross- 
talk between the amplifiers produce a series of artefacts arranged 
horizontally wi th respect to each star in the field at regular 64 
pixel intervals l lFinger et alj|2008l) . These artefacts appear crater- 
like and although are produced for every object in the field, they 
are only detectable above the noise if the star is brighter than ap- 
proximately J = 14.5 mag (although this is strongly dependent on 
the seeing). The artefacts were most pronounced in the deep J im- 
ages. 2" x 1" rectangular regions were masked at 64 pixels inter- 
vals from each star brighter than J ~ 14.5 across the entire detector 
quadrant. Cross-talk was greatly reduced in programs observed af- 
ter ESO semester 82. 
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Table 1. Summary of the target fields and the HAWK-I observations. 500 MHz radio luminosities, given in column 3, are from lMilev & De Breuck (200 



Radio galaxy 


Radio galaxy 
coordinates (J2000) 


Radio luminosity 
Log(500 MHz) 


Centre of field 
(J2000) 


Redshift 


Filter 


Exposure time 
(hours) 


5o~ depth 

(Vega) 


Seeing 
(FWHM) 


MRC 1324-262 


13:26:54.6-26:31:42 


28.46 


13:26:51 -26:30:55 


2.28 


J 


3.42 


24.9 


0.7 












H 


0.88 


23.6 


0.7 












Ks 


0.42 


22.0 


0.42 


USS 1425-148 


14:28:41.7 -15:02:28.4 


28.66 


14:28:53 -15:01:43 


2.35 


J 


3.15 


25.0 


0.65 












H 


0.72 


23.0 


0.6 












Ks 


0.56 


22.0 


0.6 


MRC 0406 -244 


04:08:51.5-24:18:16.7 


29.03 


04:08:43 -24:17:35 


2.43 


.1 


3.38 


24.8 


0.4 












H 


0.84 


23.5 


0.7 












Ks 


0.53 


22.8 


0.74 


MG 2308 + 0336 


23:08:25.1 +03:37:03.9 


28.51 


23:08:28 +03:36:15 


2.46 


J 


3.40 


25.5 


0.55 












H 


0.88 


23.4 


0.55 












Ks 


0.56 


22.4 


0.5 


MRC 2 104 -242 


21:06:58.3-24:05:09.0 


28.84 


21:06:55 -24:04:31 


2.49 


J 


3.38 


25.3 


0.65 












H 


0.67 


23.8 


0.75 












Ks 


1.53 


22.6 


0.6 


MRC 2 139 -292 


21:42:16.6-28:58:38.2 


28.74 


21:42:10-28:59:32 


2.55 


.1 


3.58 


25.3 


0.45 












H 


1.70 


24.0 


0.55 












Ks 


0.53 


22.7 


0.6 


Control field 






01:28:04-26:26:08 




.1 


3.18 


25.2 


0.7 












H 


0.84 


23.5 


0.65 












Ks 


0.53 


22.7 


0.7 



Bright stars and nearby galaxies were also masked because 
they cover a significant amount of area in some fields and there- 
fore reduce the total area where faint galaxies could be detected. 
All objects brighter than Ks = 16 mag were masked out. Finally all 
pixels with less than 30 percent of the maximum effective exposure 
time were masked as these regions are incomplete even at bright 
magnitudes and could lead to spurious detections. 



increases by this correction by a few percent, up to 20% for the 
smallest objects. 

Uncertainties on the fluxes in each band were estimated 
through a 1 = 5c74grd + °Poi S son> where °bkgrd is the image back- 
ground noise and o"p ; sson is the Poisson uncertainty in the num- 
ber of detected electrons (flux in ADU x effective instrument gain) 
measured by SExtractor. 



2.4 Source detection 

Sources were detected using SExtractor dBertin & Arnoutslll996j) 
in double image mode using weighted Ks images as detection im- 
ages. The weighted images, created by multiplying the reduced Ks 
images by the square root of their associated effective exposure- 
time maps, take into account the varying background noise level. 
Sources were defined as groups of 5 pixels that are 1.25o" above 
the local background noise. The data were filtered, before object 
detection, with Gaussian filters of width and size that matched the 
seeing. 

The colour of objects were measured from the PSF-matched 
images using apertures of approximately 2 arcsec diameterQ Total 
Ks fluxes were measured from the high-resolution Ks images using 
Kron apertures (FLUX_AUTO parameter of SExtractor). 

Photometry in Kron apertures can underestimate the total flux, 
especially if the source is faint. Therefore aperture corrections were 
applied to convert the Ks Kron fluxes to total Ks fluxes (K, ot ). Cor- 
rection factors were determined from stellar growth curves of bright 
unsaturated stars in the Ks images: a circular aperture, of the same 
total area as the Kron aperture of each object, was used to mea- 
sure the amount of light missed by the Kron aperture. The amount 
of aperture correction depends on the PSF of the image, and on 
the source size and brightness, but the source flux was typically 

1 The exception is the MRC 2 1 39 - 292 field, for which 2.2" apertures were 
used, since the stellar growth curves of the 3 images only converge to less 
than 3% beyond this diameter. 



2.5 Image completeness 

The point-source and galaxy completion of the Ks HAWK-I data 
was calculated by simulating 3000 point sources or galaxies on 
each image (in batches of 1000) using the IRAF packages GAL- 
LIST and MKOBJECT. The simulated galaxies were set to have a 
2.5 pixel half-light radius, as observed for bright galaxies in the 
images, and we assume the elliptical-to-spiral fraction is 0.4, and 
maximum elliptical axial ratio is 0.3. 

The image completion is uneven due to the varying effective 
exposure time, so each simulated source was assigned an effective 
exposure time calculated from the weights map. The sources were 
then detected with SExtractor, grouped into 4 bins of effective ex- 
posure time, and the completion level calculated for each bin. Re- 
gions with low effective exposure time have lower completion lev- 
els than regions with high effective exposure times. The average 
completion for the entire image was created by averaging the com- 
pletion of the 4 exposure-time categories, weighted by the fraction 
of image area that lies within each exposure time category. 

The average galaxy completion of each image is shown in 
Fig-H The 90% completeness is the limit to which the entire im- 
age is complete to at least 90%, including regions with low effec- 
tive exposure time. Below this limit the galaxy completeness is not 
uniform across the images because of the varying effective expo- 
sure time. The images are complete to ~ 0.6 mag fainter for point 
sources. 
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Figure 1. The average completeness of extended sources in all 7 observed 
Ks fields. The images are complete to ~ 0.6 mag fainter for point sources 
than extended sources. The entire image is complete to at least the 90% lim- 
iting magnitude, including regions with low effective exposure time. Below 
this limit the galaxy completeness is not uniform across the images. 90% 
galaxy completeness limits are given in the legend. 



2.6 Stellar classification 

At bright magnitudes (K < 19 mag), SExtractor can adequately se- 
lect unresolved objects based on their morphology, but this selec- 
tion method is not adequate at faint magnitudes. Therefore objects 
were classified as stars based on their near-infrared colours. The 
near-infrared colours of bright stars were examined in the colour- 
magnitude diagrams of the 7 observed fields, and the following 
colour-selection criteria were designed to separate stars from galax- 
ies. Stars were defined as objects which obey all of the following 
conditions: 



J-H<0.55 n H-K <0.85 n J-K<1.15 



(1) 



These criteria were tested using the Pickles stellar library. All stars 
are selected by this criteria, except for very red stars such as L- or 
early T-type dwarfs, but the density of such stars is very low. 

We identify stars in both the control fields and the HzRG fields 
using the same colour criteria so the results presented in this work 
are not sensitive to moderate changes in these criteria. Changing 
these colour cuts by up to 0. 1 mag does not affect our conclusions. 



3 IDENTIFYING OVERDENSITIES IN THE RADIO 
GALAXY FIELDS 

3.1 Galaxy number counts 

Differential if -band galaxy number counts were measured from the 
HzRG fields and the control field (CF) and are compared to galaxy 
counts from the literature in Fig.[2] The galaxy number counts were 
not corrected for completeness, and no attempt was made to correct 
for the difference in the filter bandpasses. 

The observed number counts are fully consistent with litera- 
ture counts up to the completion limit of the images: there is no 
large galaxy excess in any HzRG field compared to the control 
field or the literature. Any protocluster structure associated with 
the HzRGs may only contribute a small fraction of the galaxies in 
the observed field. 

The wide-field HAWK-I images contain a large contribution 
of foreground and background objects which dominate over any 
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Figure 2. Galaxy number counts in K for the 6 radio galaxy fields compared 
to the control field (CF) and values taken from the literature. No complete- 
ness correction was applied to the radio galaxy fields or CF. There is good 
agreement between the radio galaxy fields and data from the literature, in 
particular between 17.5 and 20.5 mag. The radio galaxy fields do not con- 
tain a large excess of galaxies. 



structure associated with the radio galaxy. Additionally, cosmic 
variance of these galaxies may hide a protocluster, or make a re- 
gion appear overdense through chance superposition. It is therefore 
necessary to remove as many interloping galaxies as possible be- 
fore searching for structures around the HzRGs. 



3.2 Selecting galaxies at 2.2 < z < 2.7 

To identify galaxy overdensities near the HzRGs we attempt to iso- 
late galaxies at 2.2 < z < 2.7. For galaxies in this redshift range 
the Balmer and 4000A break lies between the / and H filter band- 
passes. 



3.2.1 Colours of galaxies with photometric redshifts in the range 
2-2 < Zphot < 2.7 

Fig-IDplots the near-IR colours of all galaxies within t he 0.8 deg 2 
^-selected Ultra Deep Survey (UDS; lLawrencell2007t) . This field 
is covered by deep B, V, R, i', z', J, H, K, Spitzer 3.6/xm and 
4.5jUm data. An upda ted version of the A'-selected catalogue of 
IWilliamsetall j2009j) is used, in which the data from UKIDSS 
Data release 1 is su pplemented with //-ba nd data from the UKIDSS 
Data release 3 (see lWilliams et al .l2009bl for details). In addition to 
the photometric ca talogues we use the ph otometric redshift cata- 
logues compiled in|Williams et alj J2009bl) . We refer the reader to 
IWilliams et al. |{2009ab) for further details on both catalogues. 

A sample of galaxies were selected from the UDS which have 
photometric redshifts in the range 2.2 < z < 2.7. These galaxies 
have similar redshifts as the HzRGs, and are likely to have simi- 
lar colours to the companions of the radio galaxies. To obtain the 
correct galaxy colours, it is essential that the Balmer/4000A break 
of the UDS galaxies falls between the J and H bands, just like the 
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Figure 3. Near-IR colours of galaxies in the 0.8 deg 2 UDS (grey circles), 
UDS galaxies with 2.2 < z p hot < 2.7 (black solid circles), and the 6 ob- 
served HzRGs (open circles). Red lines mark the JHK criteria designed 
to select galaxies at the same redshift as the HzRGs. 98% of UDS galax- 
ies with 2.2 < z D ho t < 2.7 satisfy the ALL-JHK criterion (solid black line). 
iBruzual & Charlotl 120031) models of galaxies at z = 2.4 are overplotted. 
Each track corresponds to galaxies with ages evolving from 0.01 to 2.6 Gyr 
(J — H reddens with ages). The 6 tracks have different star formation histo- 
ries: continuous star formation (solid line), and exponentially declining star 
formation (dotted lines) with T = 1,0.5,0.1 and 0.05 Gyr (from bottom to 
top) and T = 0.01 Gyr (dashed line). The reddening vector is shown by the 
black arrow. 

HzRGs. Therefore we only select galaxies whose redshift fits have 
log y^ 1 < 2.9, and a narrow redshift probability distribution, such 
that the 3cr confidence intervals of the redshift estimate lie within 
the redshift range 2.2 < z < 2.7. This resulted in a sample of 89 
galaxies with K„„ < 20.6 mag. These galaxies are plotted in Fig.|3] 
as filled black circles. Hereafter this sample is referred to as z p hor 
UDS. 

3.2.2 Colour criteria to select 2.2 < z < 2.7 galaxies 

iKaiisawa et aflfcOOd) and lKodama et al. I42007|) devised a 2-colour 
selection criterion, 

J-H>H-K+0.5 n J-K> 1.5 {JHK). 

which selects galaxies whose Balmer or 4000A breaks are brack- 
eted by the J and H passbands, so it is well suited to our purpose 
of selecting galaxies at the same redshift as the z ~ 2.4 HzRGs. 
Galaxies selected by this cr iterion are referred to as JHK galaxies 
following the convention of Kajisawa et al. (2006). 

Approximately half (53%) of UDS galaxies who satisfy the 
JHK criteria have photometric redshifts in the range 2.2 < z p hot < 
2. 7. This selection techn ique has also been spectroscopically tested 
bv lDohertvet all ( |2010|) . who found that 10 out of 18 JHK galaxies 
lie in the redshift range 2.2 < z < 2.7. Therefore approximately 
half of the JHK galaxies in the HzRG fields are likely to lie near 
the redshift of the HzRG. 

This JHK criterion is biased toward selecting galaxies with 
low star formation rates, and cannot detect z ~ 2.4 galaxies with 
continuous star formation histories (see Fig.[3}. Over 85% of the 
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Figure 4. Comparison of the 3 control fields: the solid line is the 
230 arcmin 2 control field (CF+) consisting of the observed control field and 
the outer regions of the radio galaxy fields; the dotted line is the 56 arcmin 2 
observed control field (CF); and dashed line is 0.8 deg 2 UDS. Error bars 
come from Poisson statistics. There is good agreement between all 3 con- 
trol fields for both colour-selected galaxy populations. 

UDS galaxies with photometric redshifts in the range 2.2 < z p hot < 
2.7 do not meet the JHK criteria, so many protocluster galaxy can- 
didates may be missed. We therefore define a second colour crite- 
rion 

J-H>H-K-0A5, {ALL-JHK), 

which is marked by the solid black line in Fig.[3] 98% of UDS 
galaxies with photometric redshifts in the range 2.2 < z p hat < 2.7 
fall within this colour criterion. However only a minority of this 
sample will lie at the same redshift as the HzRGs, since only 4% of 
the UDS galaxies that satisfy this criterion have photometric red- 
shifts in the range 2.2 < Zphot < 2.7. 

Lumin ous radio galaxies are among th e brightest galaxies at 
all redshifts l lRocca-Volmerange et al.l2004l) so accompanying pro- 
tocluster galaxies are likely to be fainter than the HzRGs. The 
brightest HzRG has K tot = 17.7 mag, so we assume that all galax- 
ies with K tot < 17.5 mag are interlopers. Only the MRC2104-242 
field contains an excess of galaxies with K tot < 17.5 mag (see 
Fig.|2j. Visual inspection of these galaxies reveal that they have 
large angular sizes and blue near-infrared colours, so they are likely 
to be low redshift interlopers. For the rest of this work galaxies 
with K tot < 17.5 mag are removed from the catalogues of the HzRG 
fields and control fields. 

3.3 Comparison control fields 

Reference control fields are required to identify overdensities in 
the radio galaxy fields. A control field (CF) was observed with the 
same instrumental set-up as the radio galaxy fields, allowing for 
a direct comparison without observational bias. This control field 
is a single HAWK-I field-of-view (56 arcmin 2 ) so is sensitive to 
cosmic variance. Therefore a larger control field was compiled by 
supplementing the CF with the outer regions of the 6 radio galaxy 
fields i.e. the regions beyond 3 arcmin (~ 5 Mpc co-moving at z = 
2.4) of the HzRGs. This combined control field results in a total 
area of 230 arcmin 2 , hereafter referred to as CF+. 

The cumulative surface density within the control field (CF), 
CF+ , and the 0.8 deg 2 UDS are compared in Fig.|4] The surface 
density of galaxies in the CF+ is generally lower than or at the same 
level as the CF for all magnitudes, so the CF+ field is not contam- 
inated by large numbers of protocluster galaxies. The 230 arcmin 2 
CF+ is therefore used as the control field in the rest of this work. 

A larger control field is required to determine the uncertain- 
ties due to the field-to-field variation of galaxy number counts. 
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Figure 5. Surface overdensity within 1 .5' of the radio galaxies (top) and be- 
yond 3' (bottom). JHK galaxies are plotted in green while ALL-JHK galax- 
ies are plotted in black. The surface overdensity was measured using galax- 
ies with K lot < 20.6 and the HzRGs are not included. Three out of 6 fields 
contain significant overdensities in these two colour-selected populations. 
None of the fields are overdense beyond 3', which tells us that the overden- 
sities are not due to zeropoint errors. 



For this purpose we use the 0.8 deg 2 UDS. Colour transforma- 
tions are required to shift the observed UDS colours to HAWK-I 
colours since the UDS was observed with filters that have different 
passbands than HAWK-I filters. These transformation^ were de- 
termined using stellar population models of high-redshift galaxies, 
and by matching the galaxy number counts and colour distributions 
of the JHK and ALL-JHK populations in the UDS to the CF+ (i.e., 
ensuring the UDS cumulative number counts in Fig.|4]match those 
of CF+). Therefore these transformations may not be suitable for 
sources with different colours, such as stars. 



3.4 Surface density of galaxies in the HzRG fields 

The galaxy surface overdensity £ g is the excess surface density of 
galaxies. It is calculated as £„ = (£ fo — £) /£, where £ £ s is the ob- 
served surface density, and £ is the expected surface density mea- 
sured from the CF+ control field. 

The surface overdensity of the JHK and ALL-JHK galaxies 
was measured within 3 co-moving Mpc (1.8') of the HzRGs. This 
distance corresponds to approximately the virial radius of a mas- 
sive local cluster. The fields were selected to include the HzRGs, 
so these galaxies were not included in either the JHK or ALL-JHK 
sample. The results are shown in Fig.[5]and the number and over- 
density of each population are listed in Table [2] The error bars rep- 
resent lrj field-to-field variations in the surface density measured 
from 10,000 randomly positioned cells within the UDS. The fields 



Khawk-i = K UDS + 0.09(fl - K) UDS - 0.0635; 
(H - K) HAWK -, = 1 M(H - K) UDS - 0. 138; 
(J - H)hawk-i = 0.91(7 - H) UDS + 0.0673. 



containing USS 1425 - 148, MG 2308 + 0336 and MRC 2104- 242 
are at least Iff overdense in both JHK and ALL-JHK galaxies. The 
other 3 fields show no significant overdensity in either population. 

The surface overdensity of the JHK and ALL-JHK galaxies 
were also measured beyond 3' of the HzRGs (see bottom panel of 
Fig.|5}. There are no significant overdensities in the outer regions of 
any field so the overdensities near the radio galaxies are not caused 
by zeropoint errors or inadequate subtraction of stars. 

MRC 0406 -244 has a slight overdensity in JHK galaxies 
within 1.5' and ale under-density in the outer regions. Therefore 
the 1.8' cell around the radio galaxy is significantly overdense in 
comparison to its local surroundings suggesting MRC 0406 - 244 
may have several nearby companions. 

The spatial distribution of JHK and ALL-JHK galaxies in the 
6 radio galaxy fields are shown in Fig.(6] A visual inspection of 
these maps confirms that the fields containing USS 1425 - 148, 
MG 2308 + 0336 and MRC2104-242 contain significantly more 
ALL-JHK and JHK galaxies than expected, whilst the 3 fields con- 
taining MRC 1 324 - 262, MRC 0406 - 244 and MRC 2 1 39 - 292 do 
not. 



3.5 Significance of the surface overdensities 

The significance of the galaxy overdensity is the joint probability 
of finding an overdensity in both JHK and ALL-JHK populations. 
Since the populations are not mutually exclusive, the significance 
of each population cannot be simply combined. Instead 10,000 ran- 
dom 1.8' radius cells within the 0.8 deg 2 UDS were used to de- 
termine the probability of finding a cell which is as dense as the 
regions surrounding the HzRGs in both ALL-JHK and JHK popu- 
lations. 

The probability of finding a 1.8' radius cell within the UDS, 
with surface overdensities (£ / M ) equal or greater than each of 
the radio galaxy fields is listed in the bottom row of Table [2] 
Between 26% and 43% of all UDS regions were as dense as 
the 1.8' radius cells around MRC 1324-262, MRC0406-244 or 
MRC 2139-292, so these fields do not contain more galaxies than 
expected. The probability of finding regions as dense as those sur- 
rounding USS 1425 -148, MG 2308 + 0336 and MRC 2104 -242 
is ~ 0.3%, so the number of galaxies in these fields deviate from 
the expected galaxy density by 3o". 



4 AUTO-CORRELATION FUNCTION ANALYSIS 

If the galaxies responsible for the overdensities in the 3 overdense 
HzRG fields are physically associated they should be strongly clus- 
tered. Whereas if the overdensity is caused by a chance line-of- 
sight alignment the clustering signal is not expected to be stronger 
than average. This clustering signal was measured using the auto- 
correlation function. 



4.1 Angular correlation function 

The two-point angular correlation function, w(8), is defined as the 
excess probability of finding two objects, in regions SQ.\ and SO.2, 
separated by a distance 6, 



5P = £ 2 fa [l+w(e)]5fi 1 5fi 2 , 



(2) 



where £ £ s is the surface density of galaxies in the field consid- 
ered. Since £ j s is taken into account, w(9) will not be stronger in 
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Field 


MRC 1324-262 


USS 1425-148 


MRC 0406 -244 


MG 2308 ±0336 


MRC 2104-242 


MRC2139-292 


Type of galaxy 














ALL-JHK 


-0.2±0.2 (65) 


0.3 ±0.2 (91) 


-0.2 ±0.2 (66) 


0.4±0.2 (110) 


0.5 ±0.2 (115) 


-0.3 ±0.2 (58) 


JHK 


-0.1 ±0.7 (2) 


1.9 ±0.8 (6) 


0.3 ±0.7 (3) 


1.6 ±0.7 (6) 


1.7±0.7 (6) 


-0.2±0.7 (2) 


Fraction of 1.8'radius UDS 
cells as dense as HzRG field 


38% 


0.3% 


26% 


0.4% 


0.3% 


43% 



Table 2. Galaxy overdensities of all colour-selected populations within 1.8' of the 6 HzRGs. The number of galaxies in each sample are given in the parentheses. 
We do not include the HzRGs in these calculations. The probability of finding a 1.8' radius cell in the UDS which is at least as dense in both colour-selected 
populations is given in the bottom row. The overdensities around USS 1425 - 148, MG 2308 ± 0336 and MRC 2104-242 deviate from the expected galaxy 
density by ~ 3o~. 



MRC1324-262 MRC0406-244 MRC2139-292 
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Figure 6. Spatial distribution of the JHK (red) and ALL-JHK (black) galaxies in the 6 radio galaxy fields. Galaxies with K lot < 20.6 are marked as filled 
circles, while JHK galaxies with 20.6 < K lo , < 22.0 are open circles. The yellow star marks the positions of the HzRGs and the dotted circles denotes 1.8' 
from the HzRGs. The fields in the top row have no significant overdensity within 1.8' of the radio galaxies, whilst the fields in the bottom row contain > lo" 
overdensities of both ALL-JHK and JHK galaxies. 



overdense regions just because there are more galaxies. We use the 
lLandv & Szalavl dl993h estimator for w{8): 

DD(6)-2DR(6)+RR(e) 

w(e) = rW) ' () 

where DD(6), DR(6), and RR(d) are the number of data-data, 
data-random, and random-random pairs respectively, with angular 
separation between 6, and 6j + \. 

To effectively sample the geometry of the detection area and 
take into account boundary effects, the random catalogues were 
constructed in the following way. The weight map image was di- 
vided into 4 regions of approximately equal effective exposure 
time. The number of objects with K 101 < 20.6 mag were counted 



in each region, and the random catalogue constructed by placing 
100 times as many objects at random locations in each of the 8 re- 
gions. Terms involving random points were normalised such that 
I. e DR(8) = T. B RR(0) = I. e DD{6). 

The angular correlation function was measured in 4 indepen- 
dent bins: 0.001 - 0.1', 0.1-0.3', 0.3-0.5', an d 0.5-0.7'. The un cer- 
tainty was estimated by bootstrap resampling fcing etaljl986h and 
is quoted at the lo" level. To measure the clustering signal w(6) was 
approximated by a power law, 

w(G) = A w e~P -IC, (4) 

where A w is the amplitude of the angular correlation function, and 
p is the slope. We corrected for the integral constraint, which re- 
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Figure 7. The angular auto-correlation beyond 3' (open circles and dotted 
line) and within 1.8' (filled circles and solid line) of the 3 HzRGS in over- 
dense fields. The angular correlation function of all galaxies near the HzRGs 
is consistent with further afield (top panel). The bottom panel shows the an- 
gular correlation function of a sample of galaxies that contains 1.5 times 
more galaxies near the HzRGs than in the outer regions. These galaxies are 
more strongly clustered close to the HzRGs, where there is a significant 
overdensity, than further afield, suggesting the galaxies responsible for the 
overdensity are physically associated. 



suits from the finite field size, using a constant, IC, det ermined 
from t he ra ndom catalogues fo llowing the approach used bv llnfantej 
d 19941) and lRoche et"aTT ( l 19991) : 



IC: 



2 _ Z,A lr 0,. p RR(6i) 



(5) 



LiRR.{6i) 

The slope was assumed to be f5 = 0.6 dAdelberger et alj|2005t) as 
there is insufficient data to fit both the slope and the amplitude. The 
amplitude of the angular correlation function, A w , was determined 
by fitting Eqn.[4] 



We measured the amplitude of the angular correlation func- 
tion of galaxies with J — H>H — K + 0.25 colours within 1.8' 
of USS 1425 - 148, MRC 2104 - 242 and MG 2308 + 0336 (A imer ). 
This was compared to the A w of galaxies with the same colours and 
magnitudes that lie beyond 3' of the same HzRGs (A ou ter)- In this 
sample the clustering signal of the excess galaxies should not be 
strongly masked. If the excess galaxies are physically associated, 
we expect Aj nnel - to be greater than A oll t e r- 



4.3 Results 

The auto-correlation functions are plotted in Fig.|7]for all galax- 
ies (top panel) and galaxies with J — H > H — K + 0.25 (bottom 
panel). The auto-correlation function of the galaxies beyond 3' of 
the HzRGs and within 1.8' are plotted as open and solid circles 
respectively, and the amplitudes A oute i- and Aj nnel - are given in the 
legend. 

As expected, A w of all galaxies near the 3 overdense HzRGs is 
consistent with that measured in the outer regions, Aouter — A[ nne! - — 
0.1. The angular correlation function of both regions agree at all 
scales, suggesting there is no bias due to the different regions, and 
that w(6) and A M , for the overdense sample may be robustly com- 
pared without fear of systematic bias. 

Galaxies with y - # > # - A: + 0. 25 colours that lie within 1 . 8' 
of the HzRGs have larger w(8) at all scales compared to the outer 
region of the fields. A; nner is 2 times greater than A ou t e r, however the 
uncertainties are large due to the small number of galaxies in this 
galaxy sample, and this result is significant to only la. Nonethe- 
less the galaxies responsible for the overdensity are more strongly 
clustered near the HzRGs than similarly selected galaxies at larger 
radii. This suggests the galaxies responsible for the overdensity are 
physically associated with one another. 

We have checked our results assuming /3 = 0.8, which is 
also commonly used in the literature. We find that A w is reduced 
by a factor of approximately 2, but the disparity between A outer 
and A; nnel - is still present for galaxies with J — H > H — K + 0.25 
colours. 



4.2 Galaxy sample 

To test the angular correlation function analysis we measure 
A„. for all galaxies with 17.5 < K tot < 20.6, both within 1.8' 
(A inner ), and beyond 3' of USS 1425 -148, MRC 2 104 -242 and 
MG 2308 + 0336 (A outa ). Since the galaxies responsible for the 
overdensity make up a negligible fraction of this galaxy sample, 
any strong clustering signal of the excess galaxies will be masked, 
and we expect A; nner to be consistent with A outer . 

To recover the clustering signal of the galaxies responsible for 
the overdensity we select a sample of galaxies which contains a 
large excess. All 3 overdense fields contain ~ 3 times the expected 
density of JHK galaxies, but there are too few galaxies in this sam- 
ple to obtain a measure of w(0) in 4 independent 9 bins. We there- 
fore select galaxies that meet a similar, but slightly more relaxed 
colour criterion (J — H>H — K + 0.25). 58 galaxies with such 
colours lie within 1.8' of the 3 overdense HzRGs (a number density 
of 2arcmin~ 2 ), compared to 279 in the CF+ (1.2arcmin~ 2 ), indi- 
cating an overdensity of 0.7. Therefore the excess galaxies make up 
40% of this sample. 



5 PROPERTIES OF THE PROTOCLUSTER GALAXY 
CANDIDATES 

In this section we aim to measure the properties of all protocluster 
galaxy candidates in the 3 overdense HzRG fields. 

Only a small fraction of protocluster galaxy candidates will 
be in the JHK sample, as this criterion selects less than 15% of 
galaxies at 2.2 < z p kot < 2.7. Thus the JHK galaxies are not repre- 
sentative of the entire galaxy population surrounding the HzRGs. 

Most of the z ~ 2.4 protocluster galaxy candidates are ex- 
pected to be ALL-JHK galaxies, as this sample encompasses 98% 
of galaxies with 2.2 < z p hoi < 2.7. However the overdensity of 
ALL-JHK galaxies in the 3 overdense HzRG fields is ~ 0.4, so 
more than 70% of ALL-JHK galaxies lie in the foreground or 
background. This interloper population must be removed from the 
ALL-JHK sample to reveal the properties of the protocluster galaxy 
candidates galaxies. 
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5.1 Deriving a flux-limited sample of protocluster galaxies 
through statistical subtraction 

We do not know which particular ALL-JHK galaxy is an interloper, 
but we can statistically subtract the interloping galaxies from the 
ALL-JHK sample. 

To derived the statistical galaxy excess, the number density, 
magnitudes and colours of the interlopers were measured from the 
control field (CF+), and subtracted from the galaxy counts of the 
3 overdense HzRG fields. If the galaxy overdensities are due to 
protoclusters, this method should ideally remove all foreground and 
background galaxies, to leave only the protocluster galaxies. 

A major advantage of this method is that, by construction, a 
flux-limited sample of candidate protocluster galaxies is selected, 
irrespective of colour or star formation history. This allows us to 
determine the colours and magnitudes of all protocluster members. 

Accurate statistical subtraction requires both the control field 
(CF+) and the overdense fields be complete to the same magni- 
tude. Thus the galaxy catalogues are limited to the 90% complete- 
ness limit of the lowest-exposure time region in the shallowest field 
(K t ot < 20.6 mag), and the derived sample of excess galaxies are 
complete to this flux limit. 

5.1.1 Uncertainties and biases 

The greatest uncertainty in the statistical subtraction method comes 
from cosmic variance of the interloping galaxy population. To 
quantify the uncertainty, we ran 1000 Monte Carlo simulations us- 
ing the UDS, by selecting and combining 3 random 2.1' radius re- 
gions and subtracting the interloping galaxies as measured from a 
randomly positioned 230 arcmin 2 region. The standard deviation of 
the resulting statistical excess (or deficit) gives the 1 o" uncertainty 
in the removal of the interloping population. When deriving the 
uncertainty on a galaxy property, such as colour, galaxies were di- 
vided in colour bins and the uncertainty is the standard deviation 
within each colour bin. 

To check that selecting ALL-JHK galaxies does not bias the re- 
sults, we tested the statistical subtraction method using all galaxies 
with K to , < 20.6. The derived properties of the excess galaxies are 
fully consistent with the following presented results, which sug- 
gests that few, if any, protocluster galaxy candidates are removed 
by making the ALL-JHK selection. 

5.2 Radial extent and volume of the overdensities 

To determine the radial extent of the overdensity, we measure 
the statistical excess of ALL-JHK galaxies within a projected ra- 
dius, R, of the HzRGs. This was achieved by counting ALL-JHK 
galaxies within a radius R of USS 1425-148, MRC2104-242 
and MG 2308 + 0336, and subtracting the expected number of 
ALL-JHK galaxies within this area (determined from CF+). The 
radial extent of the overdensity is marked where the cumulative 
distribution stops increasing. 

By inspection of the cumulative distribution in Fig.[8] we find 
that the average overdensity extends to a projected radius of ~3.5' 
(6 co-moving Mpc at z = 2.4) and contains approximately 75±20 
galaxies with K, o! < 20.6. Half of the statistical excess lies within 
1.8' (3 co-moving Mpc) of the HzRGs, so this central region is 
twice as dense as the entire overdensity. 

Assuming the overdensity lies in a cube of side 12 Mpc, the 
co-moving of the overdensity is ~ 1700 Mpc 3 . Dark matter simula- 
tions of the early Universe show that matter collapses in filaments, 



co-moving Mpc 
2 4 6 




R (arcmin) 

Figure 8. The cumulative number of excess galaxies with K to t < 20.6 within 
a projected distance R of the 3 HzRGs in overdense fields (USS 1425 - 148, 
MRC2104-242 and MG 2308 + 0336). The radial extent of the average 
overdensity is located where the cumulative distribution stops increasing. 
Approximately 220 excess galaxies are situated within 3.5' (6 co-moving 
Mpc) of the 3 radio galaxies. The average protocluster contains 1/3 of the 
shown numbers. The grey shaded area shows Iff field-to-field variations. 

and our selection method preferentially selects filaments that are 
close to being perpendicular to the plane of the sky. Therefore this 
estimate of the volume is very uncertain. 

5.3 Colours of the excess galaxies 

The colour distribution of the statistical excess was derived by sub- 
tracting the colour distribution of the interloping ALL-JHK galaxies 
(measured from the CF+ and normalised by area) from the colour 
distribution of ALL-JHK galaxies within the 3 HzRGs in overdense 
fields. 

Ideally the colours of all excess galaxies within the HzRG 
fields should be derived. However, the variance due to the interloper 
population is approximately proportional to the field area. There- 
fore the colours of the statistical excess was only measured within 
3.5 co-moving Mpc (2.1') of USS 1425 - 148, MG 2308 + 0336 and 
MRC 2 104 - 242, where the number of excess galaxies is more than 
4 times the uncertainty (see Fig.[8j. 

5.3.1 Comparing the colours of the excess galaxies to control 
field galaxies 

Fig.|9] displays the colour distribution of ALL-JHK galaxies in the 
CF+, galaxies within 2.1' of the 3 overdense HzRGs, and the de- 
rived statistical galaxy excess. The colour distribution of the statis- 
tical excess differs from the control field, which is quantified using 
a two-sided Kolmogorov-Smirnov (KS) test. 

The KS test requires continuous unbinned data, so a mock cat- 
alogue of galaxy colours was created by sampling random values 
whilst constraining the catalogue to match the observed colour dis- 
tribution of the statistical excess. The KS test was performed 1000 
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Figure 10. The observed near-infrared colours and K, ot fluxes of the statistical galaxy excess (solid black histograms) compared to those of the z p h t-UDS 
sample (dashed red histogram). z p i 10l -UDS galaxies are selected from the UDS to have (i) 2.2 < Zphot < 2-7, (ii) log % 2 < 2.9, and (iii) the 3o~ confidence 
intervals of the redshift estimate lie within the redshift range 2.2 < z. < 2.7. There are 87 galaxies in the z p i, 0l -UDS sample, and the distributions have been 
normalised to the number of excess galaxies. The colours and masses of z p hoi-UDS and excess galaxies are similar with KS probability of p=0.82. The red 
tail seen in both J — H colour distributions is unique to galaxies at 2.2 < z < 2.7, whose redshifted Balmer break shifts between the J and H passbands. The 
J — H colour distribution is bimodal, with a strong peak of blue galaxies and an extended tail of red galaxies. There is no prominent peak of red galaxies. The 
expected positions of the z = 2.4 red sequence are marked as dot-dashed, dotted, and dashed lines assuming the galaxies formed at z/ = 3,4,5 respectively. 
Most protocluster galaxies are blue so they are still forming stars at z, ~ 2.4, and the red galaxies are consistent with a formation redshift of z/ ~ 3 
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J-H 

Figure 9. The top panel shows the J — H colour distribution of galax- 
ies in the 3 overdense fields (dashed histogram) and CF+ galaxies, nor- 
malised to the same area (grey histogram). The difference between these 
two distributions reveal the J — H colours of the excess galaxies, which 
is shown in the bottom panel in the solid black histogram. A two-sided 
Kolmogorov-Smirnov test comparing the colours of the excess galaxies to 
the CF+ galaxies results in prob= 0.001, implying the colour distribution 
of the two groups of galaxies differ significantly. The excess galaxies are 
therefore not randomly drawn from all redshifts, but probably comprise a 
group of galaxies at a particular redshift. 



times and in each test the colours of the excess galaxies were drawn 
from a different set of random values. The probability, p, gives the 
significance level of the test (ranging between and 1) and is the 
median of all derived p. 

The KS test rejects the null hypothesis with p=0.001, so the 
colour distribution of the excess galaxies differs significantly from 
the colour distribution of galaxies in the CF+. The excess galaxies 
in the 3 overdense fields are not drawn randomly from the control 
field population and are therefore unlikely to be caused by a chance 
line-of-sight alignment. 



5.3.2 Comparing the colours of the excess galaxies to Zphot ~ 2.4 
galaxies 

The colour distribution of UDS galaxies with photometric redshifts 
in the range 2.2 < z p i w t < 2.7 (z p hot~ UDS sample; see Section 
I3.2.U is shown in Fig.[l0] The distribu t ion is similar to that of 
2 < Zphot < 3.5 galaxies in IWuvts et al .1 d2007l) which span simi- 
lar fluxe|3- Both distributions peak at U — V ~ 0.5, and there are 
relatively few galaxies r edder than U — V = 1.0, or bluer than 
U-V = 0.0 (see Fig. 8 in lWuvts et alj2007h. The z P w-UDS sam - 
ple is bluer than the 2 < z < 3 galaxies of lvan Dokkum et "aH d2006h . 
but their study is limited to galaxies with stellar masses greater than 
10 u M o which are generally brighter than our flux-limited sample. 
We conclude that the colour distribution of the z p h ot -UDS sample 
is a good representation of z ~ 2.4 galaxy colours. 

The J — H and H — K colours, and fluxes of galaxies compris- 
ing the statistical excess are compared to the UDS galaxies with 
2.2 < Zphot < 2.7 in Fig.[To] The flux distributions of both samples 
are consistent, so a robust comparison of the two samples can be 
made. 



3 (U-V) rcs ,(AB) = 1.26 x (J-H) -0.558 for galaxies at z ~ 2.4 
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The colour distributions of the excess galaxies are very simi- 
lar to the UDS Zphot ~ 2.4 galaxies. The red tail in the J — H colour 
distribution is unique to galaxies at 2.2 < z < 2.7. It is caused by 
Balmer break of galaxies at this redshift lying between the J and 
H passbands. Both the excess galaxies and the z v hot ~ 2.4 galax- 
ies have this red tail in the J — H colour distribution. A KS test 
does not reject the null hypothesis that the excess galaxies and the 
UDS 

Zphot ~ 2.4 galaxies are drawn from the same underlying dis- 
tribution (p=0.82). Since the excess galaxies have similar/ — H and 
H — K colours, and magnitudes to z ~ 2.4 galaxies, it is likely that 
they are situated at the same redshift as the HzRGs. 



6 DISCUSSION 

The 3 HzRGs USS 1425 -148, MG 2308 + 0336 and 
MRC2104-242, are surrounded by galaxy surface overden- 
sities that deviate from the expected density by 3a. The excess 
galaxies responsible for the overdensities are strongly clustered, 
suggesting they are physically associated. One of the strongest 
pieces of evidence that the candidate protocluster galaxies are 
associated with the radio galaxies is the similarity between the 
colours of the candidate protocluster galaxies and z p hot ~ 2.4 
galaxies in the K~ selected Ultra-Deep Survey (see Fig.llOt. It 
would be a strange coincidence that the colour and luminosity 
distributions match so well if both these samples did not consist 
of z ~ 2.4 galaxies, especially considering they were selected by 
different methods. In this section we examine the properties of 
these overdensities. 



6.1 Evolution of the overdensities around z ~ 2.4 HzRGs 

6.1.1 Average mass enclosed in the overdensities 

Assuming all the mass within the volume collapses into a single 
virialized structure, the mass of the collapsed structure is estimated 
by 

M = pV(l + 5,„) (6) 

dSteidel et al.lll998h where V is the volume of the overdensity, p is 
the average density of the Universe (1.9 x 10 MoMpc -3 ), and <5„, 
is the mass overdensity, related to the galaxy overdensity through 
the bias parameter 5„, = S g /b. 

To determine the overdensity of galaxies, we compare the 
numb er density o f galax ies in the overdensities to the z = 2.5 
field. IWuvts et al.1 J2009h measure the density of z ~ 2.4 galax- 
ies with M> 4 x 1O 1O M to be 3.3 x 10~ 4 Mpc~ 3 . Stellar masses 
of the excess galaxies w ere estimated from K, ot usi ng the same 
iBruzual & Charlotl d2003l) stellar synthesis models as IWuvts et al.1 
J2009h . Six stellar synthesis models (plotted in Fig.[3j were red- 
shifted to z = 2.35, 2.46 and 2.49 (i.e., the redshifts of the radio 
galaxies). Moderate amounts of dust extinction (up to Ay=\ mag) 
were applied to the models, and the models were convolved with 
the HAWK-I J, H and Ks filter bandwidths to obtain the follow- 
ing relation: log M„ = [K tot -43.32 (±0. 15) - 2.79 (±0. 16) x (J - 
H)]/ — 2.5. The uncertainties are derived from the 3o" scatter about 
the linear relation. 

Some galaxies responsible for the overdensity have J — H ~ 
1.5 so the flux-limited sample is only complete to a mass-limit of 
~ 6 x 10 10 M G . Therefore the number density of M> 4 x 10 10 M Q 
galaxies in the HzRG fields, and the mass estimates of the overden- 
sities, are lower limits. 



The number density of excess galaxies with M> 4 x 10 Mq 
within 6Mpc of the 3 overdense HzRGs is 24 ± 12/VMpc -3 . 
The uncertainties quoted include uncertainties in the mass esti- 
mation and the statistical subtraction of the interlopers. Assuming 
V ~ 1700 Mpc 3 , the galaxy overdensity is S g ~ 42±21. The bias 
parameter for gal axies with a numbe r density of 3.3 x 10~ 4 Mpc~ 3 
at z ~ 2.4 is ~ 3 l lOuadri et al.ll2007h . thus the mass overdensity is 
S m = 14 ±7. 

The mass contained within the average radio galaxy overden- 
sity is ~ 2-7 x 1O 14 M0 which is comparable to the mass of 
present-day clusters of galaxies, and is consisten t with the mass 
estim ates of other HzRG-selected protoclusters dVenemans et al.1 
l2007h . 

The volume of the overdensity is very uncertain (see Section 
15.21 ), but this does not affect the estimation of mass enclosed within 
the overdensity, as the volume in Equation[6] cancels with the vol- 
ume used to calculate S m . 

6.1.2 Collapse time 

In a matter dominated Universe, the time for a uniform sphere to 
separate from the Hubble flow and collapse is ~ (l/6^Gp) 1//2 
jPeacocklll999h . The collapse time for these cluster-sized overden- 
sities is ~ 6 Gyr. This estimate is very uncertain as the volume of 
the overdensity is poorly constrained, but this short collapse time 
implies these overdensities may virialize by the present day. 

In conclusion, the 3 HzRG fields with significant projected 
surface overdensities fulfil the main physical requirements for 
evolving into a galaxy cluster or group by the present day. 

6.2 Rest-frame U — V colours of galaxies in z > 2 
protoclusters 

The rest-frame U — V colours of the flux-limited sample of proto- 
cluster candidates are derived through (U — V) res t(AB) = 1.26 x 
(J — //observed) — 0.558 and shown in Fig.llOl The colour distri- 
bution is bimodal, consisting of a dominant blue sequence, and a 
poorly populated red tail. 77 ± 10% of the protocluster galaxies are 
blue with (U — V)ab < 0.85. However, this is an upper limit as 
blue galaxies are detected to lower masses than red galaxies in the 
flux-limited sample. 

The blue protocluster candidates have a typical colour of U — 
Vab ~ 0.6 mag, which is bluer than local spiral galaxies, and typical 
of a present-day irregular galaxy. Such blue colours suggest that 
these galaxies are forming stars. 

The colours of the galaxies in the red tail of the distribution 
are consistent with that of a cluster red sequence which formed the 
majority of its stars at z/ ~ 3, a mere 0.5 Gyr before the epoch of 
observation. However, the red sequence is poorly populated com- 
pared to those observed in clusters at z < 1.6 (e .g.,[ Blakeslee et al.1 
120031 : IPapovich et a"L1l2010l : iTanaka et alj|20ich . The errors on the 
observed colours and the slightly different redshifts of the HzRGs 
may increase the scatter of the galaxies on the red sequence. 

A colour bimodality, and prominent red sequence has been 
found up to z ~ 2.5 using large area phot ometric surveys (e.g., 
I Williams et al. I l2009ak iBrammer et al]|2009h . These surveys probe 
much larger volumes than the protoclusters, and detect greater 
numbers of massive galaxies that populate the red sequence. The 
lack of a strong red sequence in the 3 protocluster candidates may 
be due to the low number of massive galaxies within each proto- 
cluster. 
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No significant difference is observed between the colours of 
the protocluster galaxies and field galaxies at z p hot ~ 2-4. This sug- 
gests that global environmental effects have yet to become impor- 
tant in clusters by z ~ 2.4. However, the samples of protocluster 
and field galaxies are selected by different methods, so we refrain 
fro m making deta i led co mparisons between these two populations. 

ISteidel et al.l d2005l) found a significant fraction of the galax- 
ies in a z = 2.3 protocluster were detected in the rest-frame UV. 
Hence these galaxies were still forming stars and were relatively 
blue at this redshift, which is consistent with our findings. Through 
detailed SED modelling, the ages and stellar masses of the z = 2.3 
protocluster galaxies were found to be a factor of 2 greater than 
field galaxies. While our data does not allow us to perform a sim- 
ilar analysis, future follow-up observations may reveal how wide- 
spread such evolutionary trends are among samples of protoclus- 
ters. 

Our results show that the majority of galaxies in the z > 2 
protocluster candidates have colours typical of star-forming galax- 
ies. The colour, low scatter and non-evolving slope of the tight red 
sequence in lower-redshift clusters imply the stars that constitute 
cluster ellipticals formed contemporaneou sly at z > 2 dEllis et al] 
ll997l ; lGladders et al.ll99§l ; lMei et alj2009h . The scatter and colour 
of the red sequence in a z = 1.62 cluster indicate s this epoch of 
inten se stellar-buildup occurred at z = 2.4 ±0.15 jPapovich et al .1 
whilst the evolution in the mass-to-light ratio of ellipti- 
cal galaxies up to z = 1.3 imply a formation redshift of 2 — 2.5 
dvan Dokkum & van der Mare ll l2007t) . The z > 2 protocluster can- 
didates may have been caught in this starburst phase, during which 
they form the bulk of the stellar population that will later reside in 
their elliptical galaxies. 



7 CONCLUSIONS 

We have conducted a near-infrared survey of the environment of 
6 high redshift radio galaxies (HzRG) at 2.2 < z < 2.7. We use 
colour-cuts to identify galaxies in this redshift interval, and find ex- 
cess numbers of these galaxies around 3 HzRGs: USS 1425- 148, 
MRC 2104 - 242 and MG 2308 + 0336. 

Within 1.8' (3 co-moving Mpc) of these 3 HzRGs approxi- 
mately 3 times more colour-selected galaxies are observed than ex- 
pected. The surface density of galaxies with colours of 2.2 < z < 
2.6 galaxies deviate by 3d from the field density. The clustering of 
the excess galaxies responsible for the overdensity was measured 
through the angular auto-correlation function. The excess galaxies 
are more strongly clustered than similarly-selected galaxies in the 
control field, suggesting they are physically associated. 

We derive a flux-limited sample of protocluster galaxy candi- 
dates by statistically subtracting the fore- and background sources 
from the 3 overdense fields. The observed colour distribution of 
the statistical galaxy excess is consistent with that of galaxies with 
photometric redshifts in the interval 2.2 < z< 2.7, suggesting that 
the excess galaxies are located in a narrow redshift range near the 
HzRGs, and are therefore likely to be protocluster galaxies associ- 
ated with the HzRGs. However, spectra of the protocluster galaxy 
candidates are required to confirm whether the overdensities are 
coherent structures. 

The large-scale structure around the HzRGs extends approx- 
imately 12 Mpc (co-moving), and contains 75 ± 20 galaxies with 
K tot < 20.6 mag. The average mass contained within these struc- 
tures exceeds 1O 14 M,0, and these dense regions may collapse and 



virialize by z ~ 0. Thus these objects may evolve into groups or 
clusters of galaxies by the present day. 

We find the rest-frame U — V colour distribution of the flux- 
limited sample of protocluster galaxies is bimodal, consisting of 
a dominant blue sequence and a poorly populated red sequence. 
77±10% of protocluster galaxies have blue colours, typical of star- 
forming galaxies. This implies that the bulk of stars that constitute 
cluster ellipticals have not formed by z ~ 2.4, or are too blue to lie 
on the red sequence. The U — V colours of the protocluster member 
candidates are consistent with a formation redshift of z < 3. 
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